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SUMMARY 


Table  9  summarizes  the  equation  of  state  and  mixed-ph£^  regime  parameters  for  kaolinite 
clay  and  the  pelitic  rocks  examined  in  this  study.  Examination  of  die  data  in  light  of  the 
compositions  of  the  rocks  suggests  that  the  low  pressure  properties  of  pelitic  rocks  are 
dominated  by  the  relative  fractions  of  clays  and  micas,  as  well  as  Ae  presence  of  zeolitic 
water  in  open  crystal  structures.  The  effect  of  the  clay/mica  ratio  is  primarUy  on  the 
density,  as  clays  are  less  dense  than  micas,  but  have  bulk  moduli  very  similar  to  those  of 
micas.  Zeolitic  water  has  a  significant  effect  on  both  the  density  and  the  bulk  modulus  of 
the  rocks  in  their  low  pressure  phases. 


The  high  pressure  portions  of  the  Hugoniots  represent  states  in  which  the  minerals  in  the 
rocks  have  transformed  to  high  pressure  polymorphs.  Models  of  the  high  pressure 
behavior,  obtained  by  considering  the  mixing  of  the  high  pressure  phases  of  the  individual 
minerals,  results  in  Hugoniot  states  that  are  more  dense  than  the  experimentally  measured 
states.  This  suggests  that  the  current  equation  of  state  of  one  or  more  phases  common  to  all 
three  rock  types  is  in  error.  The  high  pressure  phase  density  seems  to  be  correlated  with 
the  H2O  content  of  the  rock,  with  high  density  rocks  having  lower  water  contents, 
properties  are  affected  mostly  by  the  abundance  of  minerals  containing  no  volatiles,  relative 
to  clays  and  carbonates. 


The  use  of  a  constant  specific  heat  in  the  equation  of  state  models  allows  easy  calculation  of 
the  release  adiabats.  Comparison  with  the  release  adiabat  data  shows  that  both  kaolinite 
and  calcareous  shale  high  pressure  phases  convert  to  a  lower  density  material  upon  release 
from  the  shock  state.  Slate,  however,  may  retain  the  high  pressure  phase  to  very  low 
pressures.  The  hysteresis  in  the  shock-release  path,  and  thus  the  fraction  of  shock  energy 
deposited,  can  be  calculated  using  these  release  adiabats.  Because  of  formation  of  low- 
density  material  upon  release,  the  calculated  energy  deposition  will  be  greater  than  the 
actual  value.  Incorporation  of  such  information  into  models,  along  with  suitable  models 
for  the  low-pressure  behavior,  will  allow  the  decay  of  a  shock  wave  resulting  from 
irreversible  heating  to  be  calculated. 
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SECTION  1 
INTRODUCTION 


Clays,  shales  and  slates  make  up  a  series  of  genetically  related  sedimentary  and 
metamorphic  rocks,  called  pelitic  rocks,  characterized  by  high  contents  of  hydrous 
aluminosilicate  minerals.  Pelitic  rocks  form  an  abundant  class  of  materials  found  at  and 
near  the  earth's  surface  in  a  wide  range  of  geologic  environments.  Essentially  all  ancient 
sedimentary  basins  now  contain  shales  or  slates  and  many  modem  basins  contain 
unconsolidated  clay  mixtures  which  will,  over  geologic  time,  become  shales  and  perhaps 
slates.  Because  pelitic  rocks  are  so  widespread,  making  up  ~30%  of  all  sediments  and 
meta-sediments,  accurate  knowledge  of  their  behavior  under  shock  compression  and 
adiabatic  release  is  essential  for  the  modelling  of  shock-producing  events,  such  as  impacts 
and  explosions,  in  sedimentary  and  metamorphic  terranes. 

In  this  study,  we  conduct  an  experimental  investigation  of  the  shock  and  release  behavior 
of  a  clay  (kaolinite),  a  shale,  and  a  slate,  and  combine  our  experimental  results  with 
previously  reported  data  for  similar  materials  to  develop  equations  of  state  (EOS)  of  clays 
and  pelitic  rocks.  We  also  address  the  ability  to  predict  the  equation  of  state  of  pelitic 
rocks  from  a  knowledge  of  the  composition.  The  equation  of  state  of  kaolinite  is  expected 
to  be  very  similar  in  most  respects  to  those  of  other  clay  minerals.  The  kaolinite  EOS  is 
compared  to  that  of  muscovite  mica  to  examine  the  consequences  of  conversion  of  clays  to 
micas  during  metamorphism  of  shales  into  slates. 

We  also  use  the  equations  of  state  and  experimental  data  to  study  the  shock-release  path 
hysteresis  and  consequent  energy  deposition  in  pelitic  materials.  Knowledge  of  such 
behavior  allows  us  to  draw  conclusions  about  the  decay  rate  of  shock  waves  expanding 
from  an  impact  or  confined  explosion. 
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SECTION  2 

GENETIC  RELATIONSfflPS  BETWEEN  CLAYS,  SHALES,  AND  SLATES 


Clays  are  part  of  a  class  of  hydrous  aluminosilicate  minerals  called  phyUosilicates, 
characterized  by  sheets  of  silica  and  aluminum  hydroxide  separated  by  layers  of  cations. 
Clays,  admixed  with  other  minerals  (notably  quartz,  SiOp,  form  the  starting  material  for 

formation  of  shales  and  slates.  Shales  are  sedimentary  rocks  formed  by  the  physical 
deposition  of  clay  minerals  with  very  fme  grains  of  other  minerals  and  subsequent 
compaction  and  cementation.  The  cementing  material  most  often  is  orgarac  material  formed 
from  detrital  remains  of  organisms  or  carbonate  minerals  (usually  calcite,  CaCOg)  deposited 

by  chemical  precipitation.  As  such,  the  properties  of  shales  can  be  modelled  as  a  mixture 
of  clay,  other  detrital  minerals,  and  the  relevant  cementing  material.  One  complication  in 
the  terminology  is  that  "clay"  is  both  a  mineralogical  and  a  textural  term,  denoting  either  a 
member  of  the  clay  group  of  phyllosilicate  minerals  or  any  sediment  that  is  very  fine¬ 
grained.  The  primary  complication  is  that  the  prominent  layering  in  shales  can  introduce 
anisotropy  to  the  wave  propagation  behavior. 

Slates  result  from  relatively  low-grade  metamorphism  of  shales.  The  most  common 
metamorphic  reaction  is  the  conversion  of  clays  to  micas  with  the  evolution  of  water. 
Continued  metamorphism  to  higher  grades  results  in  the  production  of  phyllites,  which  are 
very  similar  to  slates  but  have  experienced  textural  coarsening  due  to  recrystallization  of 
minerals,  schists,  in  which  grain  growth  has  become  extreme  for  micas,  and  gneisses, 
which  represent  extreme  metamorphism,  just  short  of  wholesale  melting  of  the  rock. 
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SECTION  3 

EXPERIMENTAL  TECHNIQUES 

3 . 1  SAMPLES  AND  PREPARATION. 

Pure  kaolinite  (no  detectable  contaminants  in  X-ray  diffraction  analysis)  from  Georgia  was 
obtained  from  a  major  rock  and  mineral  sample  supplier.  These  specimens  were  compact 
but  unconsolidated  and  contained  negligible  contaminants.  The  material  was  disaggregated 
in  a  mortar  and  pestle,  moistened  with  deionized  water,  and  pressed  in  a  special  die  under  a 
pressure  of  ~1  kbar  into  pellets  37  mm  in  diameter  and  4-5  mm  thick.  The  f)ellets  had 
densities  in  the  range  of  1.93-2. 16  Mg/m^. 

A  calcareous  shale  was  collected,  by  us,  for  this  study.  It  was  obtained  from  the  lower 
Chambless  Limestone  (Cambrian)  in  the  southernmost  Marble  Mountains  of  the  central 
Mojave  Desert  in  California.  Volatile-free  oxide  composition  was  determined  via 
inductively-coupled  plasma  mass  spectroscopy  and  X-ray  fluorescence  spectrometry.  CO^ 

in  carbonates  was  determined  by  reaction  with  acid  and  volumetric  measurement  of  gas 
evolved.  The  composition  is  given  in  Table  3-1.  The  samples  had  a  mean  nonporous 
(crystal)  density  of  2.615  Mg/m^.  Ultrasonic  measurements  of  longitudinal  and  transverse 
wave  velocities  give  Vp  =  4.505  km/s  and  Vj  =  2.753  km/s,  respectively. 

Slate  samples  were  obtained  from  Portland,  Maine,  and  the  vicinity  of  Slatedale, 
Pennsylvania.  Both  slates  have  similar  properties,  with  an  average  bulk  density  of  2.756 
Mg/m3.  The  Pennsylvania  slate  has  been  studied  previously  [Furnish,  1992;  Davies  and 
Smith,  1994]  and  its  mineralogical  composition  is  reported  by  Davies  and  Smith  [1994], 
They  also  give  ultrasonic  elastic  wave  velocities  for  this  slate  as  =  4.08  km/s  and  = 
2.19  kra/s. 

3.2  EXPERIMENTAL  PROCEDURE. 

The  experimental  techniques  for  this  study  are  discussed  in  detail  elsewhere  [Ahrens, 

1987].  The  surfaces  of  a  rectangular  sample  approximately  1.3  cm  wide,  1.5  cm  long, 
and  3-4  mm  thick,  were  ground  flat  and  parallel.  In  contrast,  the  kaolinite  samples  were 
first  reduced  to  a  powder  and  cold  pressed  in  a  die.  Harness  and  parallelism  was  achieved 
to  within  ±0.02  mm.  The  sample  was  then  mounted  on  a  metal  driver  plate  (Figure  1)  with 
standard  buffers  and/or  inclined  mirrors  placed  on  the  free  (rear)  surface  of  the  sample. 
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Shock  arrival  mirrors  were  placed  against  the  rear  surfaces  of  the  driver  plate,  sample,  and 
buffers.  During  the  experiment,  the  driver  plate  was  impacted  by  a  gun-launched  projectile 
containing  a  metal  impactor  (flyer  plate).  The  impact  generated  a  planar  shock  wave  which 
propagates  through  the  driver  plate,  sample,  and  buffer  materials.  Arrival  of  the  shock 
wave  at  the  free  surface  of  each  component 
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Table  3-1.  Oxide  composition  analyses  of  calcareous  shale  from  Marble  Mountains, 
central  Mojave  Desert,  California. 


Oxide 

Weight  %  aCP)^ 

Weight  %  (XRF)b 

SiOz 

44.59 

44.20 

TiOz 

0.84 

0.76 

AlzOs 

9.68 

9.42 

FezOs 

4.77 

4.73 

CaO 

19.09 

19.30 

MgO 

1.33 

1.39 

NazO 

0.33 

0.13 

KzO 

2.04 

2.25 

COz*^ 

14.34 

14.34 

HzOd 

4.21 

3.46 

Total 

100.94 

99.98 

(a)  Inductively-couple  plasma  source  emission  spectrometry 

(b)  X-ray  fluorescence  spectrometry 

(c)  Determined  by  gas  evolution  on  treatment  with  acid. 

(d)  Determined  by  thermal  desorption  and  decomposition  (different  samples). 
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caused  disruption  of  the  specular  reflection  from  the  shock  arrival  mirrors.  One  spatial 
dimension  of  the  target,  including  all  shock  arrival  mirrors,  was  imaged  as  a  function  of 
time  by  a  streak  camera,  allowing  the  shock  wave  transit  times,  and  thus  the  shock  wave 
velocities  through  the  target  components  of  known  thickness,  to  be  determined.  The 
internal  energy  E,  pressure  P,  density  p,  and  particle  velocity  Up  of  the  shocked  sample  and 
buffer  are  constrained  from  the  conservation  of  mass,  energy  and  momentum  and  the 
requirement  that  pressure  within  the  two  materials  must  remain  in  contact  requires  that 
particle  velocity  be  continuous  across  material  interfaces. 

The  buffers  and  inclined  mirrors  are  used  to  measure  points  on  the  release  adiabats  of  the 
samples  by  the  requirement  that  the  P-u^  condition  in  the  released  state  of  the  sample  be 
identical  with  that  of  the  shocked  buffer  or  the  free  surface  in  the  case  of  the  inclined 
mirror.  Table  3-2  gives  the  properties  of  the  materials  used  for  flyer  and  driver  plates  and 
buffers  in  this  study. 

The  shock  Hugoniot  of  a  material,  which  is  the  locus  of  final  states  achieved  by  shock 
compression,  can  be  conveniently  described  as  a  polynomial  representation  of  the  shock 
wave  velocity,  in  terms  of  the  particle  velocity, 

U,  =  C^+sUp  +  s'Up  +...  (3.1) 

Here  Co  is  the  intercept  for  zero  particle  velocity  and  is  d  V/d  Up .  For  simple  non¬ 
transforming  media,  Co  is  the  zero  pressure  bulk  sound  velocity  and  S  =  [(d  Kos/dP)p=o  + 
l]/4  where  Kqs  is  the  zero  pressure  isentropic  bulk  modulus.  The  polynomial  may  be 
truncated  at  first  or  second  order  (if  material  phase  are  not  apparent).  The  pressure,  P, 
density,  r,  and  internal  energy,  E,  in  the  shocked  state  are  given  by 


Ph  ~  ^0  pQO^s^p 

(3.2) 

u, 

Ph  Poo  JJ 

(3.3) 

Us  -Up 

(3.4) 
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HerCj  the  subscripts  0  and  H  refer  to  the  unshocked  and  shocked  (Hugoniot)  states, 
respectively,  and  the  subscript  00  denotes  the  possibility  that  the  initial  sample  may  be 
porous  or  may  be  in  a  different  thermodynamic  phase  from  that  in  the  final  shock  state. 
Here  V  is  equal  to  1/p,  where  the  same  subscripts  apply.  Po  is  usually  negligible  compared 
to  Pji  and  is  thus  set  to  zero,  while  Poo  is  arbitrary  and  so  we  also  set  E  oo=  0  as  the  initial 
condition. 
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Table  3-2.  Shock  Hugoniot  parameters  for  standard  materials  used  in  this  study. 


Material  po(Mg/m3)  Co(km/s)  5  5' (s/km)  Source 


Aluminum  - 
1100 

2.707 

5.386 

1.339 

1 

Aluminum  - 
2024 

2.784 

5.330 

1.34 

2 

Tungsten 

19.224 

4.029 

1.237 

1 

Tantalum 

16.650 

3.293 

1.308 

2 

Lexan 

1.193 

2.449 

4.817 

1.907 

1.498 

0.589 

1.431 

3 

4 

5 

Polystyrene 

0.0474 

0.362 

0.9165 

0.04982  6 

(1)  Marsh  [1980] 

(2)  Mitchell  and  Nellis  [1981] 

(3)  Fit  to  data  from  Marsh  [1980]  for  Up  <  2.605km  /  s. 

(4)  Fit  to  data  from  Marsh  [1980]  for  2.605km  /  s  <  Up  <  3.456km  /  s. 

(5)  Fit  to  data  from  Marsh  [1980]  for  Up  >  3.456km  /  s. 

(6)  Calculated  from  Pqq  using  expression  from  Anderson  and  Ahrens  [1993] 
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3.3  EXPERIMENTAL  RESULTS 


3.3.1  Kaolinite. 

Table  3-3  and  Figure  3-2  present  the  experimental  results  for  kaolinite.  The  tabulated 
densities  of  the  release  states  are  obtained  using  the  Reimann  integral  formalism  and 
approximating  the  release  adiabats  as  straight  line  segments  in  the  P-V  plane  [Lyzenga  and 
Ahrens,  1978].  The  U,-Up  projection  of  the  data  (Figure  3-2a)  show  a  discontinuity 
suggestive  of  a  phase  change,  although  it  is  also  possible  that  this  discontinuity  may  be  the 
result  of  variability  between  samples.  We  fit  the  three  segments  of  the  UfU^  Hugoniot  with 
Co  -  0.820  km/s  and  .y  =  2.26  for  <  1.821  km/s,  Co  =  4.589  km/s  and  5  =  0.19  for 
1.821  <Up<  2.219  km/s,  and  Co  =  1.497  km/s  and  s  =  1.59  for  >2.219  km/s.  The  data 
of  Sekine  et  al.  [1991]  for  muscovite  (Figure  3-3a)  show  a  similar  structure,  although  they 
interpret  their  data  as  showing  no  phase  transformation  under  shock  compression.  The  P- 
V  projection  of  the  kaolinite  data  (Figure  3-2b)  show  a  sudden  steepening  at  higher  shock 
pressures,  which  is  strongly  indicative  of  a  phase  change.  In  contrast  to  the  U^-Up  data,  the 
P-V  data  for  muscovite  (Figure  3-3a)  show  no  significant  structure.  As  some  clay  minerals 
have  crystal  structures  and  chemical  formulae  similar  to  those  of  muscovite,  these  data 
taken  together  suggest  that  the  presence  or  absence  of  a  shock-induced  phase  transition  in 
clay  minerals  depends  on  the  details  of  the  composition  and  structure.  Kikuchi  et  al. 

[1993]  suggested  that  kaolinite  and  smectic  clays  do  undergo  phase  transformation, 
possibly  with  the  formation  of  a  melt  or  amorphous  material,  at  shock  pressures  of  ~30 
GPa.  The  other  possibility  is  that,  in  kaolinite,  we  are  observing  a  transformation  related  to 
the  strong  heating  in  the  samples,  which  are  -20-25%  porous.  The  muscovite  samples 
studied  by  Sekine  et  al.  [1991]  were  nonporous  and  therefore  not  heated  as  strongly. 

Based  on  the  sum  of  the  evidence  available,  we  suggest  that  relatively  H20-rich  clays 
undergo  partial  devolatilization,  with  the  liberated  H2O  acting  as  a  flux  for  the  formation  of 
a  water-rich  melt,  under  the  proper  P-T  conditions.  Here  T  is  absolute  temperature  in 
Kelvins.  The  composition  of  this  melt  and  tire  P-T  regime  under  which  this  reaction  is 
favored  cannot  be  determined  from  the  presently  available  data. 

3.3.2  Shale. 

Table  3-4  and  Figure  3-4  present  our  results  for  calcareous  shale.  The  U^-Up  data  can  be 
broken  down  into  three  regions  and  fit  by  Equation  (3.1)  truncated  to  first  order.  The  fit 
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parameters  are  Cq  =  3.803i0.097  km/s  and  s  —  1.463+0.115  for  Up  <  1.12  km/s,  Cq 
4.853±0.115  km/s  and  ^  =  0.523±0.073  for  1.12  <  <  1.98  km/s,  and  Q  = 

2.40510.139  km/s  and  5  =  1.76110.054  for  Up  >  1.98  km/s.  The  high  compressibility 
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Table  3-3.  Experimental  Data  for  Kaolinite. 
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Figure  3-2.  Experimental  data  for  kaolinite:  (a)  UrUp  projection  with  empirical  fits  to  the 
data;  (b)  P-V projection  of  data  and  Hugoniot  calculated  from  U,-Up  fit. 
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Table  3-4.  Experimental  Data  for  Calcareous  Shale. 
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Table  3-4.  Experimental  data  for  calcareous  shale.  (Continued) 
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Figure  3-4.  Experimental  data  for  calcareous  shale:  (a)  U^-Up  projection  with  empirical  fits 
to  the  data;  (b)  P-V  projection  of  data  and  Hugoniot  calculated  from  U^-Up  fit 
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shown  at  low  pressures  in  the  P-V  data  and  the  definitely  defined  segments  of  the  U^-Up 
Hugoniot  suggest  that  the  lowest  pressure  data  represent  Hugoniot  states  entirely  within  the 
low-pressure  phase  of  the  shale  and  the  highest  pressure  data  represent  states  entirely 
within  a  high  pressure  phase.  The  intervening  region  represents  a  mixed-phase  regime 
where  the  formation  of  the  high  pressure  phase  is  incomplete,  so  that  the  Hugoniot  states 
fall  on  a  mixing  line  between  the  metastable  extensions  of  the  low  and  high  pressure  phase 
Hugoniot  curves.  We  should  note  here  that  the  term  "phase"  is  not  strictly  correct,  since 
the  rock  is  a  mixture  of  several  different  mineral  phases.  Rather,  the  term  "low  pressure 
phase"  in  the  present  context  is  used  to  refer  to  a  material  composed  of  the  same  mineral 
phases  as  the  original  sample,  whereas  the  term  "high  pressure  phase"  refers  to  a  material 
in  which  some  or  all  of  the  constituent  minerals  have  transformed  to  high  pressure 
polymorphs.  The  term  "mixed  phase"  in  the  present  paper  refers  to  a  material  in  which 
minerals  have  partially  transformed  to  high  pressure  polymorphs,  but  the  transformation  is 
incomplete. 

In  the  lowest  pressure  experiments,  the  plastic  shock  wave  in  the  shale  was  preceded  by  an 

elastic  precursor.  If  we  assume  that  the  velocity  Ud  of  elastic  precursor  should  be  identical 
to  Vp  under  infinitesimal  strain,  then  the  elastic  precursor  obeys  U^,  =  4.505  +  L2«^  for  U^i 

and  Up  in  kra/s. 

3.3.3  Slate. 

Table  3-5  and  Figure  3-5  give  the  present  data  for  slate.  The  two  slates  studied  are  very 
similar,  with  the  slate  from  Maine  being  slightly  more  dense.  The  experimental  data  for 
the  two  rocks  are  statistically  indistinguishable,  so  we  have  chosen  to  combine  the  two  into 
a  single  data  set.  The  same  Pennsylvania  slate  as  used  in  the  present  experiments  has  been 
studied  by  Davies  and  Smith  [1994]  and  Furnish  [1992].  To  increase  the  number  of  data 
available  and  thus  tlie  quality  of  the  results,  we  have  combined  those  experimental  results 
with  the  present  data.  They  are  also  presented  in  Figure  3-5. 

The  U  -u  data  for  slate,  as  for  the  calcareous  shale,  can  be  divided  into  three  regions 

S  p 

characterized  by  straight  line  segments.  The  parameters  for  these  are  Q  =  4.203  km/s  and 
5  =  2.837  for  Up  <  0.286  km/s,  Q  =  4.795  km/s  and  ^  =  0.768  for  0.286  <  <  2.169 

km/s,  and  Cq  =  3.220  km/s  and  ^  =  1.494  for  Up  >  2.169km/s,  for  poo  =  2.756  Mg/mk  As 
with  the  shale  data,  we  interpret  the  suaicture  in  the  slate  Hugoniot  as  indicating  the 
existence  of  single  "phase"  regions  at  high  and  low  pressures,  with  an  intervening  mixed- 
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phase  region  characterized  by  partial  transformation  of  the  low  to  the  high  pressure  phase. 
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Table  3-5.  Experimental  Data  for  Slate. 
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(.007)  (.004)  (.010)  (.070)  (.64)  (.034)  (.079)  (1.52)  (.125) 

(Lexan  Buffer) 

6.276  2.49  3.604 

(.092)  (.18)  (.125) 

(Polystyrene  Foam  Buffer) 


Table  3-5.  Experimental  data  for  slateXContinued) 
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data;  (b)  P-V projection  of  data  and  Hugoniot  calculated  from  Us-Up  fit 
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We  detected  an  elastic  precursor  to  the  plastic  shock  in  our  lowest-velocity  shot  (Figure  3- 
6),  indicating  that  material  strength  is  significant  at  low  shock  stresses.  Davies  and  Smith 
[1994]  also  state  that  their  data  show  evidence  of  an  elastic  precursor  to  the  plastic  shock 
wave  at  low  pressures.  This  suggests  that  the  low-velocity  segment  of  the  U^-Up  Hugoniot 
presented  above  is  representative  of  the  elastic  wave,  rather  than  the  plastic  wave.  Support 
for  this  conclusion  can  be  found  in  the  fact  that  the  intercept  of  the  low  velocity  segment  is 
close  to  the  measured  longitudinal  wave  velocity  of  4.08  km/s  [Davies  and  Smith,  1994]. 
If  the  Hugoniot  represented  a  purely  plastic  wave,  the  intercept  should  be  the  bulk  sound 
velocity  of  3.27  km/s. 
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Figure  3-6.  Streak  record  from  shot  963  showing  arrival  of  elastic  and  plastic  shock 
waves  at  the  free  surface  of  the  sample. 
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SECTION  4 

EQUATIONS  OF  STATE 


4.1  KAOLINITE. 


While  the  experimental  U,-Up  data  are  useful,  they  are  valid  only  for  the  initial  bulk  density 
of  the  samples  on  which  the  experiments  were  performed.  In  the  case  of  kaolinite,  this  is 
only  75-80%  of  the  theoretical  bulk  density  of  the  material.  In  order  to  extrapolate  the  data 
to  other  densities,  we  require  an  appropriate  general  thermodynamic  equation  of  state  for 
each  phase  or  mixture  of  phases.  These  are  obtained  by  fitting  the  P-V  Hugoniot  data  in 
each  single  phase  region  of  the  Hugoniot.  The  pressure  on  the  Hugoniot  can  be  related 
to  the  pressure  and  energy  on  an  isentrope  of  the  same  phase  at  the  same  density  by 


5  'y» _ 

YsizYil-Yjl 
2  Yh 


(4.1) 


where  is  the  pressure  on  the  isentrope  at  specific  volume  is  the  internal  energy 

gained  by  compression  to  volume  from  some  reference  condition  (usually  1  bar  and  298 
K)  along  the  isentrope,  y  is  the  thermodynamic  Gruneisen  parameter,  and  is  the 
internal  energy  difference  between  the  Hugoniot  phase  at  that  reference  condition  and  the 
initial  state  of  the  sample.  For  the  present  discussion,  we  will  always  use  1  bar  and  298  K 
as  the  reference  conditions  for  the  isentrope. 


There  are  a  number  of  parameterizations  that  can  be  used  for  the  reference  isentrope.  We 
have  chosen  to  use  the  3rd  order  Birch-Murnaghan  equation  of  state  (BMEOS)  because  it 
has  been  shown  to  provide  a  good  description  of  geologic  materials  and  is  widely  used. 
With  the  BMEOS,  P^  and  E^  are  given  by  [Ahrens,  1979]: 
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where 


x  = 


(4.4) 


and 


^  =  ^{A-K')  (4.5) 

where  Po  is  the  density,  Vq  =  1/po,  and  and  K'  are  the  isentropic  bulk  modulus  and  its 
pressure  derivative,  all  at  the  reference  condition.  The  Griineisen  parameter  is  usually 
expressed  in  the  fonn 


7=70 


r  y\’‘ 


v^y 


(4.6) 


For  the  present  study,  we  have  fit  the  EOS  parameters  to  the  experimental  data  using  a 
nonlinear  least  squares  fit  directly  to  Equations  4. 1-4.6.  The  value  of  Po  for  the  low 
pressure  phase  is  set  to  the  density  based  on  the  crystallographic  parameters  of  kaolinite. 
The  value  we  use  is  obtained  from  the  molar  mass  and  volume  stated  by  Robie  et  al. 

[1978].  Because  we  have  only  a  few  data,  we  have  chosen  to  assume  7o  =  1  and  n  =  1, 
since  the  fit  is  insensitive  to  the  values  of  these  parameters  and  the  assumed  values  are 
typical  for  most  materials.  The  fit  values  of  Kqs  and  K  for  the  low  pressure  phase  are  49.4 
GPa  and  2.36,  respectively. 

For  the  high  pressure  phase,  we  again  assume  7  =  1  and  n  —  We  also  assume  £„  =  0.5 
MJ/kg,  which  is  in  the  middle  of  the  range  of  values  displayed  by  the  high  pressure  phases 
of  silicate  minerals.  We  also  used  the  second  order  (in  strain)  BMEOS  (see  Ahrens 
[1987]),  in  which  ^  =  0  (i.e.,  K  =  4),  then  Ps  =  Kos  (X’  -  X5)/2.  Under  these 
assumptions,  the  fit  gives  po  =  3.359  Mg/m^  and  Kqs-  157  GPa. 


4.2  PELmC  ROCKS. 

As  much  as  possible,  we  have  tried  to  make  use  of  supplemental  data  to  constrain  the 
values  of  parameters  that  are  only  weakly  dependent  on  fits  to  the  shock  wave  data  or  when 
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the  number  of  available  data  are  limited.  We  use  measured  Archimedian  densities  for  the 
values  of  po  of  Ae  low  pressure  phases.  We  also  use  compositional  data  to  estimate  some 
properties.  By  determining  the  mineral  abundances  in  the  rocks,  we  can  use  simple  rules  to 
obtain  several  important  EOS  parameters.  The  assumption  we  make  is  that  the  process  of 
physical  mixing  without  chemical  reactions  or  phase  changes  conserves  the  total  volume  of 
material.  Under  this  assumption,  we  can  get  density,  bulk  modulus,  and  thermal 
expansion  coefficient  via 


Po 


V' 

I  y 


^0  “  Po 


Ttl; 


Poi^Oi . 


-t 


m;ao. 


i  Poi 
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(4.8) 


(4.9) 


where  m  is  mass  fraction,  a  is  volume  thermal  expansion  coefficient,  and  the  subscript  i 
refers  the  /'th  mineral  phase.  Unsubscripted  variables  refer  to  the  whole  rock.  Although 
we  could  use  values  of  po  and  Kos  from  (4.7)  and  (4.8)  as  final  EOS  parameters,  the  fits  to 
the  shock  wave  data  are  very  sensitive  to  these  parameters  and  it  is  best  to  allow  them  to 
vary  during  the  fitting  process.  Strictly  speaking,  (4.8)  only  refers  to  the  isothermal  bulk 
modulus,  but  the  difference  between  the  isothermal  and  isentropic  bulk  moduli  at  ambient 
conditions  is  negligible.  We  also  note  that  energy  is  additive,  so  that  specific  heat  and 
transition  energy  are  given  by 


I 

E„='Lm,E„  (4.11) 

i 

For  the  present  investigation,  we  have  assumed  Cv  =  3i?/p,  where  p  is  the  mean  atomic 
weight  and  R  is  the  ideal  gas  constant.  We  can  combine  the  values  of  po,  Kqs  cCq,  and  Cvo 
obtained  from  (4.7)-(4.10)  to  obtain  Yq  from 
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Note  that,  although  the  specific  heat  in  (4.12)  is  Cp^  not  Cv,  the  distinction  between  these  is 
negligible  at  ambient  conditions,  where  %  is  defined.  We  also  chose  to  estimate  n  by  an 
average  weighted  by  }f/po,  under  the  assumption  that  this  quantity  is  influenced  most  by 
phases  with  large  volumes  and  large  values  of  %.  In  all  the  EOS  fits,  we  have  held  Yo,  n, 
and  E^ai  the  values  obtained  by  the  prescriptions  described  above. 

4.2.1  Equation  of  state  of  shale. 

The  shale  we  studied  is  relatively  carbonate  rich  (32.6%  CaCOs  by  weight).  The  carbonate 
content  of  shales  can  be  quite  variable,  so  we  have  concentrated  here  on  a  carbonate-rich 
shale  and  use  published  data  for  carbonate-poor  shales  to  obtain  two  different  equations  of 
state.  We  used  the  composition  (Table  3-1),  on  an  H20-free  basis,  to  calculate  the 
equivalent  mineral  mixture  model  of  the  rock,  and  then  accounted  for  H2O  in  the  analysis 
and  corundum  in  the  noim  by  converting  feldspars  and  other  minerals  to  clays  and  micas. 
We  represent  HjO  tied-up  in  the  structure  of  clays  and  micas  as  a  separate  constituent. 

Table  4-1  gives  the  resulting  approximate  mineral  composition.  The  rock  is  composed  of 
roughly  equal  amounts  of  quartz,  clays,  and  calcite.  This  result  is  confirmed  by  powder  X- 
ray  diffraction  analysis.  The  mean  Archimedian  value  for  po  of  the  low  pressure  phase  is 
2.615  Mg/m^  with  calculated  values  of  %  =  0.73,  n  =  1.52,  Cv=  1396  Jkg-^K-k  We  held 
Kos  at  30.9  GPa,  calculated  from  (4.8),  and  find  that  the  fit  to  the  low  pressure  data  not 
affected  by  an  elastic  precursor  gives  IC  =  8.15.  The  calculated  parameters  for  the  high 
pressure  phase  are  Yo  —  1-2,  n  =  1.49,  E,r=  .42  MJ/kg,  and  Cv=  1396  Jkg'^K"'.  A  fit  to 
the  four  highest  pressure  shock  compression  data  gives  Po  =  3.659  Mg/m^,  Kos=  265  GPa 
and  K'  =  3.14. 

To  study  relatively  carbonate-poor  shales,  we  have  combined  the  data  of  Trunin  et  al. 
[1988]  and  Al'tshuler  and  Pavlovskii  [1971]  for  Russian  shales  and  of  Olinger  [1978]  for 
kerogen-free  Chemung  Shale.  Trunin  et  al.  [1988]  give  a  chemical  analysis  of  their  shale, 
which  we  use  as  representative  of  all  the  samples,  since  the  three  data  sets  are  the  same 
within  the  scatter  of  the  data.  The  analysis  give  ranges  of  abundances  of  different 
constituents,  so  that  our  estimation  of  the  mineral  content  of  the  shales  is  only  approximate. 
They  also  do  not  report  an  H2O  analysis  but,  based  on  the  normative  corundum  abundance 
and  the  densities  of  the  rocks,  the  dominant  phyllosilicates  are  probably  micas,  rather  than 
clays  with  the  rock  containing  only  -8%  CaCOa.  Table  4.2  gives  the  estimated  mineral 
composition  of  these  shales.  We  estimated  po  for  the  low  pressure  phase  on  the 
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assumption  that  all  the  nonnative  corundum  was  in  muscovite  and  get  po  =  2.813  Mg/m^. 
Other  calculated  parameters  are  %  =  0.79,  n  =  1.73,  and  Cv=  1159  Jkg-'K-'.  To  avoid 
possible  problems  caused  by  elastic  precursors  in  the  original  data,  we  fit  the  low  pressure 
phase  EOS  to  the  data  in  the  pressure  range  from  10  GPa  to  16  GPa  and  get  Kos=  61.1  GPa 
and  K  =  4.52.  Calculated  parameters  for  the  high  pressure  phase  are  Yo  =  0-93,  n  =  1.92, 
Cv=l  159  Jkg"‘K'* ,  and  =  0.62  MJ/kg.  A  fit  to  data  in  the  pressure  range  from 
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Table  4-2.  Estimated  mineral  composition  of  carbonate-poor  shale  studied 
by  Trunin  et  al.  [1988]. 


Mineral  Weight  Percent 


Quartz 

40 

Plagioclase 

15 

Pyroxene 

13 

Calcite 

8 

Muscovite 

20 

Magnetite 

4 

55  GPa  to  120  GPa  gives  po  =  3.740  Mg/m^,  Kos=  184  GPa  and  iT  =  3.24. 


4.2.2  Equation  of  state  of  slate. 

Although  we  used  slate  from  two  separate  locations,  the  total  data  set,  including  the  scatter 
in  the  data,  shows  that  the  two  are  indistinguishable  from  the  standpoint  of  their  shock 
wave  properties.  The  Pennsylvania  slate  used  in  this  study  has  also  been  studied  by 
Furnish  [1992]  and  by  Davies  and  Smith  [1994],  Table  4-3  gives  the  mineral  composition 
of  the  slate,  as  reported  by  Davies  and  Smith  [1994]  based  on  X-ray  diffraction  analysis. 
We  assume  that  the  material  reported  as  amorphous  in  the  analysis  is  actually  clay  minerals 
too  disordered  to  give  good  X-ray  diffraction  peaks.  Based  on  the  densities  and  porosities 
reported  by  Furnish  [1992],  the  value  of  Po  for  the  low  pressure  phase  is  2.823  Mg/m^ 
Based  on  the  mineral  composition,  we  get  Kqs  =  50.5  GPa,  %  =  0.71,  n  =  1.65,  and  Cv= 
1225  Jkg-iK-h  Most  of  the  low  pressure  phase  shock  compression  data  are  complicated  by 
the  presence  of  an  elastic  precursor,  seen  in  our  data  and  those  of  Davies  and  Smith  [1994]. 
Based  on  the  limited  reliable  data,  K'  =  8.0  for  the  low  pressure  phase.  For  the  high 
pressure  phase,  we  calculate  %  =  0.83,  n  =  1.8,  Cv=  1225  Jkg'^K'h  and  E,r  =  0.62 
MJ/kg.  A  fit  to  the  shock  wave  data  above  50  GPa  gives  po  =  4.110  Mg/m^,  Xo5=  432 
GPa  and  K  =  1.99. 
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Table  4-3.  Mineralogy  of  Pennsylvania  slate  reported  by  Davies  and  Smith  [1994]. 


Mineral 

Quartz 

Plagioclase 

Calcite 

Ferroan  Dolomite 
Pyrite 
Chlorite 
mile  •+■  Mica 
Amorphous 
Total 


Weight  Percent 
37 
12 
12 

7 
3 

8 

16 

5 

100 
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SECTION  5 
DISCUSSION 


5. 1  MIXED-PHASE  HUGONIOTS 

A  characteristic  of  the  shock  Hugoniot  curves  of  most  geologic  materials  seems  to  be  the 
transformation  of  the  material  to  a  high  pressure  phase  assemblage  at  shock  pressures 
typically  below  50  GPa.  A  complete  description  of  the  shock  compression  behavior  of 
these  materials  requires  that  a  suitable  description  be  found  for  the  behavior  of  this  mixed 
phase  region.  We  have  taken  the  relatively  simple  assumption  that  the  mixed-phase  region 
between  the  low  and  high  pressure  phase  regions  can  be  modeled  as  a  mixture  of  the 
Hugoniot  states  of  the  two  phases.  In  particular,  we  use  the  same  volume  conservation 
prescription  used  earlier  (Equation  (4.7))  to  estimate  po?  except  that  we  now  apply  it  to  pn 
at  a  given  value  of  P^: 


(5.1) 

where ///pp  is  the  mass  fraction  of  the  high  pressure  phase.  The  starting  state  for  both 
single  phase  Hugoniot  curves  must  be  the  same.  In  reality,  the  prescription  in  (5.1)  is  only 
approximate,  since  the  Hugoniot  temperatures  for  the  two  phases  at  the  same  single-phase 
Hugoniot  pressure  are  different,  sometimes  by  the.  order  of  10^  K.  However,  since  thermal 
expansion  (a)  is  a  relatively  small  effect,  which  become  less  significant  with  compression, 
made  smaller  by  the  decrease  in  a  as  Cf  material  is  the  effects  of  the  temperature  difference 
on  the  densities  is  much  smaller  than  the  density  difference  between  the  two  phases. 

In  addition  to  the  transformation  from  low  pressure  phases  to  high  pressure  phases, 
materials  also  eventually  undergo  melting  as  a  result  of  shock-induced  heating.  Our  data  do 
not  give  clear  evidence  for  this  phenomenon.  As  we  noted  earlier,  the  ’’high  pressure 
phase”  of  kaolinite  may  actually  be  a  mixture  of  volatile-rich  melt  and  volatile-depleted 
solid,  based  on  previously  reported  results  for  samples  which  were  shocked  and  then 
recovered.  Unfortunately,  such  studies  introduce  an  ambiguity  as  to  when  the  melting 
actually  occurs  (i.e.,  upon  compression  or  release).  Previous  work  with  silicate  minerals 
that  transform  to  high  pressure  phases,  such  as  quartz,  suggests  that  diaplectic  glasses  may 
form  when  a  high  pressure  phase  is  kinetically  inhibited  from  reverting  to  the  low  pressure 
phase  upon  release  from  the  shocked  state.  The  only  clear  evidence  of  a  further  phase 
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change,  probably  melting,  as  the  shock  pressure  increases,  is  from  the  data  of  Trunin  et  al. 
[1988]  for  shale.  Their  highest  pressure  data  (in  the  range  from  199  to  450  GPa)  show  a 
different  slope  in  the  plane  from  those  of  the  lower  pressure  data  we  have  used  in  this 
study.  This  break  in  slope  probably  indicates  that  these  data  are  in  either  a  mixed-phase 
region  or  in  new  single  phase  region  of  the  Hugoniot. 

To  model  the  mixed-phase  regions  of  the  Hugoniots  for  the  materials  considered  in  this 
study,  we  have  assumed  that  the  transformation  is  kineticaUy  inhibited,  so  that  the 
transformation  rate,  and  thus  the  fraction  of  material  transformed,  can  be  approximated  by 
an  Arrhenius  relation  of  the  form 


(5.2) 

where  £'a  is  an  activation  energy.  We  assume  that  some  phase  boundary  must  be  crossed 
before  any  transformation  takes  place.  Under  the  assumption  that  the  rate-limiting 
temperature  must  be  the  preexisting  phase,  we  use  the  Hugoniot  temperature,  Tji,  of  the 
low  pressure  phase.  is  given  by 


where  the  first  term  is  the  temperature  increase  along  the  reference  isentrope  and  the  second 
term  is  the  temperamre  increase  associated  with  the  entropy  increase  induced  by  the  shock 
wave.  As  previously  stated,  we  use  To  =  298  K  for  all  cases  considered  in  this  study. 

For  each  material  considered  in  this  study,  we  fit  the  parameters  of  (5.3)  to  the  value  of 
fiipp,  obtained  by  finding  the  value  that  gives  the  proper  density  for  each  experimental 
point,  and  the  calculated  value  of  7)/  for  the  low  pressure  phase.  We  assume  that  the 
kinetics  are  not  dependent  on  the  pressure,  only  on  the  temperature  of  the  low  pressure 
phase.  This  approach  has  obvious  shortcomings,  but  insufficient  data  are  available  to 
develop  a  more  rigorous  thermodynamic  model  with  phase  boundaries,  unless  a  fully 
thermodynamically  consistent  model,  treating  each  mineral  phase  separately  and  accounting 
for  temperature  differences  between  the  low  and  high  pressure  phases,  is  attempted.  That 
is  well  beyond  the  scope  of  this  study  which  examines  whole-rock  properties  without  the 
complications  of  considering  each  mineral  species  separately.  We  also  are  forced  at  the 
present  to  define  the  phase  boundary  between  the  low  pressure  phase  and  high  pressure 


35 


phase  as  a  constant  pressure.  While  this  must  be  in  error,  data  for  a  range  of  different 
initial  samples  densities  would  be  required  to  map  out  the  P-T  behavior  of  the  boundaries. 
Above  the  boundary  pressure,  transformation  to  the  high  pressure  phase  is  allowed,  but  not 
below. 

Based  on  the  present  EOS  results  and  the  assumption  of  the  Dulong-Petit  limit  for  the 
specific  heat,  we  get  A  =  8  and  30.5  kJ/mol  for  kaolinite,  with  the  phase  boundary  at 
12  GPa;  A  =  7.8  and  Ea  =  15.8  kJ/mol  for  the  carbonate-poor  shales,  with  the  phase 
boundary  at  17  GPa;  A  =  58  and  Ea  =  28.6  kJ/mol  for  the  calcareous  shale,  with  the  phase 
boundary  at  18  GPa;  and  A  =  1.34  and  £„  =  5.26  kJ/mol  for  slate,  with  the  phase  boundary 
at  21  GPa.  Of  course,  the  values  of  A  allow  the  value  of  fupp  to  be  greater  than  unity  above 
a  given  temperature.  Obviously,  the  value  of fijpp  must  be  set  equal  to  1  for  all  cases  when 
(19)  would  give /hpp  >  1.  Figures  5- 1-5-4  present  the  Hugoniot  curves,  obtained  using 
the  EOS  parameters  obtained  in  this  study  and  the  mixed  phase  region  description  presented 
here,  in  comparison  to  the  experimental  data. 

5.2  COMPOSITIONAL  TRENDS  IN  PROPERTIES. 

One  of  the  most  pressing  needs  when  considering  the  shock  wave  properties  of  geologic 
materials  is  to  be  able  to  use  compositional  and  metamorphic  history  information  to  estimate 
parameters  that  would  otherwise  have  required  the  performance  of  shock  wave 
compression  experiments.  Part  of  the  purpose  of  the  present  study  is  to  use  experimental 
data  to  draw  conclusions  concerning  the  effects  of  composition  and  metamorphism  on  the 
equations  of  state,  so  that  trends  could  be  identified. 

We  approach  this  problem  in  two  ways.  First,  we  can  compare  the  high  pressure  region  of 
the  Hugoniot,  calculated  from  the  Hugoniots  of  the  individual  mineral  constituents  via  the 
volume  conservation  prescription  already  presented  (i.e.,  (11)  but  mixing  Hugoniot 
volumes  at  a  given  pressure).  These  calculated  curves  are  presented  in  Figures  (5-2)-(5-4). 
The  calculated  curves  underestimate  the  volume  of  the  Hugoniot  state.  This  suggests  that 
the  high  pressure  behavior  of  some  constituent  common  to  both  shales  and  the  slate  is  at 
present  incorrectly  described.  It  may  also  indicate  that  the  amount  and  behavior  of  H2O  is 
are  being  improperly  accounted  for. 

The  second  approach  is  to  simply  try  to  correlate  trends  in  the  fit  EOS  properties  with 
trends  in  composition.  The  clearest  trend  is  in  the  high-pressure  phase  value  of  po,  which 


36 


increases  with  decreasing  H2O  content.  Although  slate  apparently  does  not  follow  the  trend 
established  by  the  other  materials  studied,  we  must  remember  that  the  lack  of  HjO  analyses 
for  the  carbonate-poor  shale  and  the  slate  make  estimates  at  to  water  content  highly 
uncertain. 

One  of  the  most  important  conclusions  to  be  drawn  concerns  the  H2O  content  of  the  rock. 
Politic  rocks  in  which  the  dominant  phyllosilicate  minerals  are  micas  (including  chlorite  for 
this  discussion),  tend  to  be  more  dense  than  those  containing  the  more  clays  such  as 
kaolinite,  which  contain  more  structural  water.  This  can  be  traced  directly  to  the 
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Calculated  Hugoniot  and  release  adiabats  for  calcareous  shale,  compared  with  the  experimental 


Figure  5-3.  Calculated  Hugoniot  and  release  adiabats  for  carbonate-poor  shale,  compared  with  the 
experimental  data. 
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Figure  5-4.  Calculated  Hugoniot  and  release  adiabats  for  slate,  compared  with  the  experimental  data. 


principle  of  volume  conservation  under  mechanical  mixing  and  the  greater  density  of  the 
less  water-rich  phyllosilicates.  The  presence  of  zeolitic  water,  as  was  found  in  the 
calcareous  shale  studied  here,  has  an  important  effect,  as  it  implies  the  present  of  very  open 
mineral  structures  and  thus  low  density.  These  open  structure  minerals  also  are  very 
compressible,  so  that  the  whole  rock  bulk  modulus  is  decreased  markedly,  at  least  in  the 
low  pressure  phase.  Surprisingly,  the  bulk  modulus  of  the  low  pressure  phase  of 
kaolinite,  and  presumably  other  clays,  is  very  similar  to  that  of  muscovite  (and, 
presumably,  other  micas).  Thus,  the  primary  effect  of  water  content  on  the  bulk  modulus 
is  seen  only  when  there  is  a  large  zeolitic  water  content. 

The  effects  of  the  clay/mica  ratio  on  the  densities  of  the  high  pressure  phases  is  less 
important,  since  the  clays  seem  to  undergo  a  phase  transformation,  while  the  micas 
apparently  do  not.  The  content  of  quartz  and  other  nonhydrous  minerals  has  the  most 
pronounced  effect  on  the  bulk  modulus  of  the  high  pressure  phase.  Most  of  these  minerals 
undergo  transformation  to  high  pressure  phases  with  large  densities  and  large  bulk  moduli, 
so  that  the  corresponding  rock  high  pressure  phase  will  be  relatively  dense,  as  seen  in  the 
data  for  the  carbonate-poor  shale,  which  appears  from  the  chemical  analysis  to  have 
significant  abundances  of  plagioclase  feldspar  and  pyroxene. 

5.3  RELEASE  ADIAB  ATS. 

Figures  5-1, 5-2,  and  5-4  also  present  the  experimental  P-V release  adiabat  points  obtained 
as  described  earlier.  In  all  cases,  the  release  adiabats  from  the  high  pressure  region  of  the 
Hugoniot  parallel  the  P-FHugoniot  at  high  pressures,  with  deviation  to  larger  volumes  at 
lower  pressures.  Release  from  the  mixed-phase  regions  seems  generally  to  result  in  large 
volume  increases  probably  related  to  reversion  to  the  low  pressure  phase.  Both  kaolinite 
and  calcareous  shale  show  large  volume  increases  upon  release,  especially  in  the  mixed- 
phase  region,  similar  to  the  carbonate  rocks  and  minerals-[Vizgirda  and  Ahrens,  1982; 
Anderson  and  Ahrens,  1995].  The  large  volume  increases  at  low  pressures  strongly 
suggest  that  the  high  pressure  phases  of  these  materials  revert  to  the  low  pressure  phase  or 
to  some  other  large- volume  phase  upon  release.  The  release  data  for  slate  are  somewhat 
ambiguous,  showing  either  release  paths  which  parallel  the  Hugoniot  at  high  pressures,  or 
which  show  large  volume  increases  at  high  pressures.  Because  thermodynamically 
possible  release  adiabats  should  be  straight  or  concave  upward  in  the  figure,  we  suggest 
that  the  data  showing  the  large  volume  increases  are  in  error. 
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The  experimental  release  points  can  be  used  to  check  the  validity  of  our  equations  of  state. 
The  points  shown  are  actually  lower  Uraits  on  V,  since  the  straight-line  approximation  of 
Lyzenga  and  Ahrens  [1978]  produces  the  lower  extremum  of  thermodynamically  allowable 
volume  increase  upon  adiabatic  release  with  no  phase  change.  For  a  rigorous  comparison, 
one  would  need  to  calculate  the  particle  velocity  increase  for  a  test  release  path  using  the 
Riemann  integral  and  compare  that  to  the  observed  particle  velocity  increases.  However, 
since  the  release  state  particle  velocities  often  have  fairly  large  uncertainties,  the  straight-line 
approximation  of  the  release  path  is  sufficient  for  the  present  study.  The  assumption  that 
Cv  is  constant  has  the  advantage  of  making  the  pressure  Pr  of  a  point  on  the  release  adiabat 
an  analytic  function  of  the  release  state  volume  Vr  and  the  original  shock  pressure  P a'- 

P,  =  P»  +  (Ph  -  Ps)(V„/ V,)"-'  exp[jj"  (Y/V)dv]  (5.4) 

Figures  5-1— 5-4  show  these  release  adiabats,  calculated  under  the  assumption  that  the  high 
pressure  phase  is  "frozen  in"  at  the  shock  state  mass  fraction  by  kinetics.  Reversion  of  the 
high  pressure  phase  to  some  other  lower  density  phase  is  manifested  in  the  figures  by 
experimental  release  paths  at  lower  pressures  that  are  shallower  than  the  calculated  paths. 
This  phenomenon  is  significant  for  both  kaolinite  and  shale.  However,  the  higher  pressure 
data  generally  show  good  agreement  with  the  calculated  release  adiabats. 

5 . 4  IRREVERSIBLE  HEATING  AND  ENERGY  DEPOSITION. 

One  important  use  of  calculated  release  adiabats  is  the  estimation  of  the  energy  deposited 
irreversibly  in  the  rock  due  to  hysteresis  in  the  shock-release  path.  In  terms  of  the  internal 
energy  characteristic  of  the  shock  wave,  the  fraction  f  of  the  energy  not  recovered,  and 
therefore  irreversibly  deposited  as  heat,  is  given  by 

r=i-^  (5.5) 

Eh 

where 

=  +  (5.6) 
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where  AE  indicates  the  energy  of  transformation  of  any  phase  changes,  at  the  conditions 
under  which  the  phase  changes  take  place,  and  the  subscript  ad  refers  to  an  adiabatic  path, 
which  is  isentropic  for  any  given  single  phase. 

Figure  5-5  summarizes  the  values  of  f  for  the  shales  and  slate  considered  in  this  study , 
using  the  calculated  release  adiabats  shown  in  Figures  5-2— 5-4.  Because  of  the  tendency 
for  the  real  release  adiabats  to  expand  to  larger  volumes  at  low  pressures,  these  curves  give 
erroneously  high  values  of/.  The  exception  is  for  slate,  which  seems  to  behave  as  if  the 
high  pressure  phase  persists  to  low  pressure  upon  release. 


44 


1.  8’0  9'0  fO  Z 


45 


Carbonate-Poor  Shale 


Table  5-1.  Equation  of  state  parameters  and  mixed  phase  regime  Arrhenius  equation  paramters  for  kaolinite,  shales,  and  slate. 
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